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Electric double layer supercapacitors are a fast-rising class of high-power energy storage devices 
based on porous electrodes immersed in a concentrated electrolyte or ionic liquid. As yet there is 
no microscopic theory to describe their surprisingly large capacitance per unit volume (volumetric 
capacitance) of ~ 100 F/cm^, nor is there a good understanding of the fundamental limits on volu- 
metric capacitance. In this paper we present a non-mean-field theory of the volumetric capacitance 
of a supercapacitor that captures the discrete nature of the ions and the exponential screening of 
their repulsive interaction by the electrode. We consider analytically and via Monte-Carlo simula- 
tions the case of an electrode made from a good metal and show that in this case the volumetric 
capacitance can reach the record values. We also study how the capacitance is reduced when the 
electrode is an imperfect metal characterized by some finite screening radius. Finally, we argue that 
a carbon electrode, despite its relatively large linear screening radius, can be approximated as a 
perfect metal because of its strong nonlinear screening. In this way the experimentally-measured 
capacitance values of ~ 100 F/cm"^ may be understood. 



I. INTRODUCTION 

The present energy crisis has created a growing de- 
mand for efficient, portable, and high-power energy stor- 
age devices. Electric double layer (EDL) supercapacitors 
are fast emerging as a promising potential solution to this 
problem [l|. In an EDL supercapacitor, energy is stored 
at the interface between an electron-conducting (metal- 
lic) electrode and an electrolyte or ionic liquid via the 
reversible adsorption of ions onto the electrode surface. 
In this way, counterions adsorbed onto the charged elec- 
trode effectively comprise the second half of a parallel- 
plane capacitor whose thickness is equal to the radius a/2 
of the ions. If the charge of these ions is described as a 
uniformly charged plane, as in the mean-field approach, 
then one arrives at a capacitance C which is equal to 

Ch = 2eaeA/a, (1) 

a result first envisioned by Helmholtz in 1853 0. Here, 
£0 is the vacuum permittivity, e is the dielectric constant 
of the ionic solution, and A is the total surface area of 
the electrode. In mean-field theories of the EDL, Ch/A 
plays the role of a maximum possible capacitance per 
unit area. As an example, for e = 2 and a = 1 nm Eq. 
dH) gives Ch/A ^3 ^F/cm^. 

For practical applications, a supercapacitor is best 
characterized not by its capacitance per unit area, C/A, 
but by its capacitance per unit mass or per unit volume 
("volumetric capacitance"). For this reason, there has 
been much emphasis on the development of conducting 
materials with very high specific surface area that can 
be used as electrodes. Among the more promising can- 
didates are highly porous carbons 0-0] and carbon nan- 
otube "forests" [Gj. In such devices the specific surface 
area S can be as high as § = 1000 m^/cm'^. To under- 
stand how this is possible, one can imagine an electrode 
with slit-like pores of width d = 1 nm separated by con- 
ducting walls with thickness b — 1 nm (Fig. [ij . For the 
sake of argument, we take this electrode to be the anode; 




FIG. 1. A schematic picture of the cross-section of a highly 
porous supercapacitor electrode. The solid, metallic electrode 
(gray area) has deep, planar pores of width d, separated by 
walls of thickness b. The electrode is open on one side to an 
ionic liquid (IL) or a concentrated electrolyte. 



one can imagine that the cathode is its mirror reflection 
to the right. In a supercapacitor device, the anode and 
cathode are electrically isolated by a membrane that is 
penetrable to the ions, so that the well-conducting ionic 
liquid between them forms EDLs on the tortuous surfaces 
of both electrodes. In this way the supercapacitor con- 
sists of two double layer capacitors in series; in this paper 
we are concerned with calculating the anode capacitance. 

If the electrode in Fig. [T] is placed in contact with 
an ionic solution with e = 2 and a = 1 nm, then the 
Helmholtz expression of Eq. ([T|) predicts a volumetric ca- 
pacitance C ~ 30 F/cnr^. In fact, capacitance values 
as large as C = 100 F/cm^ have been reported for such 
devices d, Q • How is this possible? 

In order to resolve this puzzle, let us briefly return to 
the problem of a planar, nonporous double layer capaci- 
tor. It has been shown recently [8] that the capacitance 
per unit area of an EDL is not necessarily limited by the 
Helmholtz value. When the charge on a planar electrode 
is small enough that adsorbed ions are separated from 
each other by a distance much larger than their diameter 
a, the mean-field approach fails and the effects of elec- 
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FIG. 2. Schematic depiction of anions (black circles) neutral- 
izing a positively charged nanopore in a metallic anode (solid 
gray area), (a) When the pore width is large compared to 
the ion diameter, opposite walls of the pore have independent 
EDLs and every anion can be said to have a single positive 
image charge (white circles with dashed outline), (b) When 
the pore width is comparable to the anion radius, anions form 
a two-dimensional (2D) charged layer within the pore. Each 
anion has an infinite series of image charges (white and dark 
gray circles with dashed outlines), which produce an expo- 
nential interaction between neighboring anions. For clarity of 
illustration, image charges are shown for one anion only. 



tronic polarization of the electrode surface must be taken 
into account. In particular, when the electrode is made 
from a good metal, each ion forms an image charge in 
the electrode surface. The ion and its image charge to- 
gether make an electric dipole which repels adjacent ions 
by a screened \/r^ interaction rather than the normal 
1/r interaction. Such a reduced interaction, along with 
the positional correlations between adsorbed ions, allows 
the capacitance of a single interface to be as much as 
three times larger than Ch in practical situations. The 
crucial importance of image forces for the structure and 
capacitance of the EDL has been recognized b y a number 
of previous authors (see, for ex amp le, Refs. MWu and the 
very recent publication of Ref. Il2h . 

The notion of a double-layer comprised of ion-image 
dipoles is also relevant for describing porous, metallic su- 
percapacitor electrodes (Fig.[T]), provided that the width 
d of the pores is much larger than the ion diameter a, 
so that opposite walls of a pore have independent, non- 
interacting EDLs [see Fig.t^Ia)]. For electrodes with such 
wide pores, enhanced capacitance can be explained using 
the theory of Ref.[^. However, in supercapacitors where d 
is comparable to a, EDLs on opposite walls of a nanopore 
merge and new physics should emerge. Indeed, recent ex- 
periments by Gogotsi and co-workers have demonstrated 
a surprising increase in the capacitance as the width of 
pores in a carbon-based electrode is made comparable to 
the diameter of bare ions in an organic electrolyte [13i | or 
in an ionic liquid Q. 

In order to explain these results, one can try to ex- 
tend the Helmholtz mean-field approach to the case of a 
narrow pore by replacing the charge of ions in the pore 
by two identical, coinciding, uniformly-charged planes lo- 
cated midway between the pore's two walls (see the result 
of a similar approach for cylindrical pores in Ref. 14). In 
this picture, each of the charged planes forms a Helmholtz 



capacitor with one of the pore's metal walls, so that the 
total capacitance is 2Chi 3.S it would be for a much wider 
pore. Thus, the mean- field approach cannot capture the 
unique effect of narrow pores. 

Going beyond the mean-field level, however, one can 
recall that charges confined within a narrow nanopore 
create an infinite series of image charges in the two con- 
ducting electrode walls [Fig. [2jb)]. This leads to an in- 
teraction between ions which decays exponentially with 
the distance between them [T5|. 

In a recent paper, Kondrat and Kornyshev [l6j recog- 
nized that such an exponentially suppressed interaction 
can lead to large capacitance C > Ch in the nanopore 
for reasons similar to those discussed in Ref. H. How- 
ever, in calculating the capacitance, the authors of Ref. 
[TgI imagined that the pore is filled with many anions and 
cations and they used a description where the small net 
charge is spread equally among all ions in the pore. The 
total electrostatic energy was then calculated using the 
exponential interaction evaluated at the average distance 
between all ions. Such an approach can be called a semi- 
mean-field approximation and does not lead to quantita- 
tively correct results, as we will show below. 

In the present paper, as in Ref. [^, we completely aban- 
don the mean-field approach in order to address a fun- 
damental question: how large can the volumetric capac- 
itance of an EDL supercapacitor be? We construct a 
theory which takes into account correlations between dis- 
crete anions and the screening of the Coulomb interac- 
tion by the conducting electrode surface. We consider 
explicitly two cases for the electrode material. First, we 
examine the case where the electrode can be considered 
a perfect metal (Sec. |TT]) and we verify our theoretical 
predictions with a simple Monte Carlo (MC) simulation 
(Sec. imp . We find that under realistic circumstances 
the capacitance per unit area can be up to ten times the 
Helmholtz value, with the corresponding volumetric ca- 
pacitance as large as 150 F/cm^. Secondly, we analyze 
the case where the electrode is an imperfect metal with a 
screening radius comparable to the pore width (Sec. llVt . 
The crossover between this theory and that of Sec. |TT] is 
carefully discussed. In Sec. |V] we consider how our the- 
ory applies to supercapacitors made with graphitic car- 
bon electrodes, and we argue that at not-too-small volt- 
ages the graphite is well approximated as a good metal 
due to the strong effect of nonlinear screening. We close 
in Sec. IVII by briefly examining a different model of a 
porous supercapacitor, where the electrode is made from 
a random assembly of conducting spheres that are three- 
dimensionally connected, and show that very large volu- 
metric capacitance can result in this situation as well. 



II. CAPACITANCE OF A SINGLE 2D METAL 
PORE 

In this section we consider an electrode made from a 
perfect metal which has deep, planar pores of width d 
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(Fig.[T]). Such pores are assumed to be in contact with an 
ionic liquid described by the restricted primitive model: 
a neutral mixture of hard-core monovalent ions with the 
same diameter a. We assume that a < d, so that ions 
can enter the pores. We also assume that d — a <C a, so 
that ions in the pore can be described as a 2D liquid. A 
voltage source provides the positive potential difference 

V between the electrode and the bulk of the ionic liquid 
that attracts anions into the pore. If Q is the amount 
of electronic charge that has moved through the voltage 
source onto the electrode relative to the state a.t V — 0, 
then the differential capacitance of the EDL is defined as 
C = dQ/dV. 

In principle, at = the pore may already contain 
some finite and equal number of anions and cations. In 
this reliable analytical calculation of the total elec- 

trostatic energy U{Q), which is necessary for calculating 
the capacitance, is very difficult. Therefore, the effect of 
allowing both ionic species to simultaneously enter the 
pore is examined only numerically at the end of Sec. IIIII 

This paper concentrates instead on the case when the 
pore is empty at = 0. This situation results when 
the chemical potential of ions in the ionic liquid is lower 
than the free energy per ion in a filled, neutral pore. 
Such a difference in chemical potential can arise from 
two sources. First, ions in a three-dimensional (3D) ionic 
liquid are surrounded by a larger number of oppositely- 
charged neighbors, which lowers the interaction part of 
the chemical potential outside the pore. Second, when 
the width d of the pore is close to the ion diameter a, the 
entropic contribution to the chemical potential inside the 
pore increases sharply. For the case of electrolyte solu- 
tions, there is also a positive contribution to the chemical 
potential associated with the necessity of stripping the 
solvation shell from each ion that enters the pore. 

We therefore assume that the pores are empty at 

V — 0. As the voltage is increased from zero, the pores 
remain empty until some finite voltage V — Vt- At 

V > Vt the pores of the anode begin to fill with an- 
ions, while cations remain away from the anode. This 
picture allows us to formulate a simple analytical calcu- 
lation of the total energy [/, presented below, based on 
the repulsion between anions in the pore. 

Our general approach to calculating the capacitance is 
as follows. We first describe the total electrostatic energy 
U (n) associated with the lowest energy configuration of 
n anions per unit area in the pore. If entropic effects are 
ignored, then the value of the charge Q of the pore is that 
which minimizes the system's total energy U—QV, where 
the term —QV represents the work done by the voltage 
source. Using the equilibrium condition (i(C/—Qy)/c?(5 — 
along with Q = eAn gives 



V 



dU 1 dU 



(2) 



dQ eA dn 

The differential capacitance of the pore C = (dV/dQ)^^ 



can therefore be written 
C = e^A 



2,2 (fU 

\ dn? 



(3) 



The capacitance can be expressed as a function of volt- 
age, C{V), by combining Eqs. © and ([3]). 

In the remainder of this section we first calculate the 
capacitance of the pore in the zero temperature limit and 
then estimate the effect of the ions' finite thermal energy. 

We begin our theoretical description by noting that a 
point charge e located in the plane halfway between two 
metal walls creates an electric potential within that plane 
equal to [13 



<t>{r) 



neoed 



J2KaH2n^l)r/d]. 



(4) 



Here, r is the radial distance from the point charge, d is 
the distance between the metal walls (the pore width), 
and Kq{x) is the zeroth order modified Bessel function 
of the second kind. At distances r > d, Eq. ^ can be 
expanded to lowest order to give 



(A(r) 



2^/2exp[—Trr/d\ e 
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(5) 



Since the sub-leading-order term of Eq. is exponen- 
tially smaller than that of Eq. ([5]), this approximation 
has a negligible effect on the capacitance and we use Eq. 
^ everywhere in further calculations. 

When a given area density n of anions is inside the 
metal pore, the repulsive interaction between anions in- 
duces strong positional correlations. In their lowest en- 
ergy configuration, the anions form a strongly-correlated 
liquid, reminiscent of a 2D Wigner crystal, where anions 
are separated from their nearest neighbors by a well- 
defined spacing ^ n^^/^. In such an arrangement the 
total repulsive energy among anions is minimized while 
maintaining the area density required to neutralize the 
electrode. 

If wc postulate a crystalline arrangement of the anions, 
then the electrostatic energy U of this state can be calcu- 
lated exactly by making use of the interaction potential 
in Eq. ([5]). Due to the short-ranged nature of the inter- 
action, this energy is well approximated by considering 
only nearest-neighbor interactions in a square lattice of 
anions. Such an approach gives 



U{n) = 2nAe(t){n ^/^) - (/i - u)nA 

= A^^{nd'f'^e^Y> 
— (/i — u)nA. 



ineosd 



(6) 



The term — (/i — u)nA takes into account the voltage- 
independent energy associated with bringing each anion 
from the bulk of the ionic liquid into the pore; fi is the 
chemical potential of ions in the bulk of the ion liquid 
and u is the self-energy of an anion in the pore. The 
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term —{fi — u)nA is linear in n and therefore, by Eq. 
([3]), disappears from the capacitance. Its only effect is 
to produce a finite threshold voltage Vt — —{^ — u)/e 
required to bring anions into the metal pore, as discussed 
above. Our theory treatment assumes that — (/i — m) > 0. 

Taking the derivative dU/ dQ as in Eq. ((2]) , we find an 
expression for the voltage in terms of the ion density: 



V^Vt 



2'K\/2 exp[~7T / VmP] 



'i-K£r,ed 



(7) 



Similarly, the capacitance can be evaluated by Eq. ([3]), 
which gives 



C 



nd 



{ndy/^eiip[n/Vnd^]CH. 



(8) 



In the limit n ^ l/d^, Eqs. Q and ([8|) can be combined 
to give an analytical expression for the capacitance as a 
function of voltage at small V — Vt: 



C ~ 327r^ 



e/Aneoed 
V -Vt 



In" 



8n' 



e/Aireoed 
V-Vt 



-,C 



H- 



(9) 

For larger voltages corresponding to {V — 
Vt) / {e / Aneoed) > 0.1 the capacitance is well-described 
by the power law relation 



,e/4:TreosdJ d 



(10) 



Eqs. ([8]) and (|9]) suggest that at low ion density (or 
small V ~ Vt) the capacitance can be much larger than 
the Helmholtz value. This result can be understood phys- 
ically by noting that at such low ion densities the frac- 
tional coverage of excess ions on the electrode surface 
na^ ^ 1, so that it is incorrect to think of the EDL in 
the mean-field way: as a uniform layer of surface charge. 
Rather, the neutralizing ionic charge consists of discrete 
ions whose interaction is exponentially small due to the 
aggressive screening by the metal pore. Positional corre- 
lations among these ions help them to avoid each other, 
resulting in a lower energy than what is possible in mean- 
field descriptions of the EDL and therefore in larger ca- 
pacitance that is not limited by the physical distance a/2 
between the electrode and its countercharge. With grow- 
ing ion density, the capacitance decreases, until at some 
finite voltage Vmax thedensity of ions in the pore reaches 
its steric limit: n ~ 1/a^. By Eq. (O, 

j2a e 

V,nax ^ 27rW — exp[-7ra/d]- -. (11) 

V a 47reo£" 

Fig. [3] shows the capacitance as a function of voltage, 
C{V), plotted for the cases d = a, 1.5a, 2a. 

If the width of pores in the electrode is increased, the 
capacitance decreases, as shown in Fig. [3] In the limit 
where the pore thickness d ^ a, as in Fig. Ufa), the ca- 
pacitance can be described using a theory of independent 



EDLs comprised of ion-image dipoles. Such an approach 
gives C « 1.3Ch per interface Q at na^ = 1 (the rela- 
tively fiat tail of the C-V curve) and Cmax ~ 5Ch per 
interface at n — >■ 0, so that the total capacitance per pore 
is smaller than the result shown in Fig. [3] by more than 
two times. As mentioned above, this "anomalous" in- 
crease in the capacitance for narrow pores is the result 
of the strong, exponential screening that results from the 
presence of two close metal walls [Fig. [2jb)]. 

Formally, Eq. ([5]) diverges as the density of ions van- 
ishes {V — Vt goes to zero). Of course, this expression 
neglects entropic effects among the ions, which are im- 
portant in the limit where ions in the pore are so sparse 
that their typical interaction energy is smaller than the 
thermal energy fc^T. At such low densities the corre- 
lated, lattice- type structure of ions in the pore disappears 
and we obtain a finite capacitance at {V — Vt) —> 0. 

In order to estimate the value of this capacitance max- 



imum, we note that when eq 



-1/2 



) < fcsT the total 



free energy F can be written using a truncated virial ex- 
pansion: 



Fc^ F,d + An^kBTB{T) - eAn{V ~ V). 



(12) 



Here, Fid = AnkBTln{na^) is the free energy of a two- 
dimensional ideal gas and B{T) is the second virial co- 
efficient. B{T) is calculated from the interaction energy 
e(/)(r) between two ions [Eq. ([5])] as 



B{T) = i 1 (^1 - exp 



(r) 



knT 



2Trrdr (13) 



167r 



(T* 



(14) 



Here, T* is defined as the dimensionless temperature 



e'^/A-KEoed 



(15) 



As an example, room temperature corresponds to T* « 
0.04 for an ionic solution with dielectric constant e = 2 in 
a pore with d = 1 nm. At temperatures that are not very 
large, T* < 0.5, the virial coefficient B{T) is larger than 
the physical area 7ra^/4 occupied by each ion, so that the 
hard-core interaction between ions is unimportant for the 
virial expansion. 

As before, we can use the equilibrium condition 
dF/dn = to give a relation between the voltage and 
the ion density n: 



e{V - Vt) = keT [2nB{T) - In (l/na^)] 



(16) 



The capacitance can also be related to n according to 
C = e'^A^{d'^F/dn'^)^^ , which gives 



C 



nd^ 



T*[l + 2nB(T)]d 



(17) 



According to Eq. ((TC)) . in the limit {V — V() = the 
ion density approaches n = Wq[2B{T) / a?]/2B{T), where 
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Wq [x] is the principle branch of the Lambert W function 
(Wo[a;] Ki \nx for x ^ 1). Over the experimentaUy rel- 
evant range of temperature 0.03 < T* < 1, the value of 
Wo[2B(r)/a2] « 1, so that Eqs. (HH) and (HH) can be 
combined to give the following approximate relation for 
the capacitance at = Vt as a function of temperature: 



Cmax (-^) 



7 



T* In^ [167r/(T*)2] d"-^' 



(18) 



Here 7 is a numerical constant; Eq. suggests 7 = 
47r^ « 39, while MC simulations (see Fig. |4] in the fol- 
lowing section) give 7 — 34 ± 1. This is a surprisingly 
good agreement, considering that Cmax is determined by 
a relatively large ion density n ~ 1/B(T), which is at 
the limit of applicability of the truncated virial expan- 
sion of Eq. (IT^ . Indeed, Eq. is applicable only in 
the "gas phase" corresponding to nB{T) <^ 1, which is 
realized at V < Vt- On the other hand, Eq. fails 
completely in the high-density correlated liquid phase, 
where nB(T) ^ 1 and the potential energy of repulsion 
between anions [see Eq. ©] dominates the entropic con- 
tribution to the free energy. A more complete theory of 
the capacitance at finite temperature and large ion den- 
sity would require a theory of the free energy of ions in 
the liquid state and cannot be captured by the truncated 
virial expansion presented here. The zero temperature 
analytical result of Eq. ([TU]), however, should be accu- 
rate in the limit where the thermal energy ksT is small 
compared to the typical interaction energy e(j){n^^^^). At 
room temperature and for d = 1 nm and e = 2, this cor- 
responds to moderately large ion density n(P > 0.3. 
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FIG. 3. (Color online). The capacitance in units of Ch of 
a 2D metal nanopore with width d, plotted as a function of 
dimensionless voltage (lower a^xis) according to Eqs. and 
(121). Different curves are labeled by their corresponding value 
of d/a and are truncated at their corresponding value of Cmax 
as given by Eq. (jlSP : here, T*a/d = 0.04. Vmax is shown for 
d/a = 1, and the top axis shows the ion concentration for this 
same case. 

The above results for the capacitance of a single metal 
pore can be used to calculate the total capacitance of the 
electrode by multiplying by the number of pores in the 



electrode. Thus, the total volumetric capacitance C is 
given by C = (eo£S/a) • C/Ch- As an example we can 
consider the electrode depicted in Fig. [U with d = 6 = 1 
nm, in contact with an ionic liquid with e = 2 and a = d. 
For such a capacitor the mostly flat tail of the C~V curve 
of Fig. 131 where C /Ch ~ 4.5, corresponds to C w 90 
F/cm'^. If the electrode can be treated as a perfect metal, 
then such a capacitor would demonstrate a peak in the 
differential capacitance at a particular voltage Vt, as in 
Fig. El For T* = 0.04, <C^ax ~ 150 F/cm^. 



III. MONTE CARLO SIMULATION OF A 2D 
METAL PORE 

In order to verify the theoretical predictions of the pre- 
vious section, we perform MC simulations of a 2D metal 
pore open to a reservoir of positive and negative hard- 
sphere ions. Our general approach is to use the Grand 
Canonical Monte Carlo (GCMC) method to impose a dif- 
ference in chemical potential between positive and neg- 
ative ions in the system, thereby simulating an applied 
voltage V . Specifically, the chemical potential of each ion 
type is specified according to 



(19) 



where /i is the chemical potential of the reservoir. The 
number of each ion species is allowed to fluctuate with 
time. We use our simulation to measure the resulting 
equilibrium number of positive and negative charges in 
the pore at a given V , which defines the net charge QiV). 
The capacitance of the pore is calculated by the discrete 
derivative dQ/dV. 

The details of our simulation method are as follows. 
We begin each simulation by randomly placing 100 of 
each type of ion on a square 2D plane of area A = 
20 X 20d^ and stipulating the dimensionless tempera- 
ture T* [see Eq. (flSl) ] and the dimensionless voltage 
V* — V / [e/ Ane{)£d). The ion diameter a is taken to 
be equal to the pore width d. Before any data is taken, 
ions are allowed to take 10^ GCMC steps to reach equi- 
librium. The number of positive and negative ions, M_|_ 
and respectively, are then averaged over the fol- 

lowing 30 to 50 X 10^ GCMC steps. The charge of the 
electrode is defined as Q = e(M_ — M_|_) and the ca- 
pacitance is given by the discrete derivative C{V) ss 
[Q[V + Ay/2) - Q{V - Ay/2)]/Ay. Care is taken to 
ensure that all results are independent of the initial ion 
configuration. 

Following the standard GCMC procedure T8], one 
GCMC step consists of either an attempted move by a 
randomly-chosen ion or the attempted addition or re- 
moval of an ion from the system. Attempted moves occur 
more often than attempted addition/removal at a ratio of 
9:1. We give the simulation area periodic boundaries, 
so that an ion leaving one edge enters at the opposite 
edge. The total electrostatic energy Utot of a given con- 
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figuration of ions is calculated as 

Utot = ^ X! *9ie0(rij), (20) 

where r^j is the distance between ions i and j (found us- 
ing the minimum image convention 19]), qi = ±1 is the 
sign of ion i, and 4)(r) is the interaction law given by Eq. 
([S])- Attempted moves and addition/removal events are 
accepted and rejected based on the corresponding change 
in Utot, as given by the traditional acceptance rules for 
GCMC [l^. The ions are treated as hard spheres, so 
that only those moves/additions/removals resulting in 
non-overlapping ions are accepted. For a more detailed 
discussion of the GCMC method see Refs. [lH andl20i. 

Fig. U shows the capacitance measured by our GCMC 
simulation as a function of voltage for a system with 
Vt* = Vt/ie/iTTEoed) = -(/x* - u*) = 0.3, calculated 
at four different dimensionless temperatures. In this sit- 
uation, the pore is essentially empty of ions at = 
and at positive voltages contains only one ionic species, 
so that its capacitance is well described by the analyti- 
cal theory of Sec. |IT1 Both the voltage and temperature 
dependence of the capacitance correspond closely to an- 
alytical predictions. Larger voltages could not be exam- 
ined by our simulation since these correspond to large 
ion fillings na^ ^ 1 (see the top axis of Fig. [3]), at which 
the simulation fails to reach equilibrium in a reasonable 
amount of time. 
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FIG. 4. (Color online). MC results for the capacitance C, 
in units of Ch, of a 2D metal pore with Vt = 0.3, plotted 
as a function of voltage for a range of temperatures T*. The 
thick line shows the prediction of our analytical theory. For 
each temperature, the capacitance attains its maximum at 
V* ~ Vt* — 0.3. These maximum capacitances (shown in 
the inset by filled squares) can be described by Eq. (|18[l with 
7 = 34 ± 1 (the theoretical curve with 7 = 4%'^ is shown by 
the thin line). 

So far we have dealt only with pores that are empty 
a,t V = 0. In the remainder of this section we use the 
GCMC simulation to qualitatively examine a pore con- 



taining a substantial, neutral concentration of both an- 
ions and cations at V = 0. To arrive at such a situation, 
one should increase the chemical potential /j,, causing the 
value of Vt to decline [2l| . 

Fig. [S] shows, as an example, the capacitance of a 
pore with fi* — u* = 0, which corresponds to moder- 
ately large ion filling ai V = Q. We also show a sys- 
tem with Vt = 0.3 for comparison. The pore with 
II* — u* =0 (triangles) is more than half-filled at zero 
vohage: (M+ -I- M_)a'^/L'^ sa 0.63. As the vohage is 
increased from zero, cations are driven out of the pore 
and anions are attracted to the pore until at V* > 0.35 
only anions remain in the pore and the capacitance is 
reasonably well described by our analytic treatment of 
the previous section (as shown by the solid line). At 
V* < 0.35, on the other hand, the strong attraction be- 
tween cations and anions affects the capacitance. As a 
rough approach to explaining this data, one may imag- 
ine that at small voltage the net ionic charge consists of a 
small number Q/e oi "excess anions" on the background 
of a large number of neutral, tightly-bound cation-anion 
pairs. These excess anions seek to maximize their dis- 
tance from each other by forming a correlated, Wigner 
crystal-like arrangement in a way that is similar to the 
description of the previous section. Under this descrip- 
tion, one may expect the same analytical theory to hold 
as for large negative fi* , since the neutral pairs are essen- 
tially non-interacting and therefore play only a small role 
in determining the capacitance. Fig. [5] suggests that this 
approach gives a reasonably accurate description of the 
finite temperature truncation of the capacitance diver- 
gence. Indeed, the capacitance in limit V* + {fi* —u*) = 
is very similar for the two C-V curves. However, this 
approach does not explain the weak capacitance maxi- 
mum at V* « 0.2 in the C-V curve corresponding to 
fi* — u* ~ 0, which remains a puzzle. 




0.1 0.2 0.3 0.4 0.5 0.6 

voltage, V/(e/4m ea) 



FIG. 5. (Color online). MC results for the capacitance of a 
2D metal pore with T* = 0.04 a.t n* — u* — (triangles) and 
— (/.I* — u*) = 0.3 (circles). The thick line is the prediction of 
Eq. ^ with V* = -(^i* - M*) = 0. 

This qualitative explanation of the filled pore data is 
similar to that of Ref. but is fundamentally different 
from the approach of Ref. [l^ These authors assumed 
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that the charge of excess anions is spread equally over all 
ions in the pore, with each ion getting a small fraction Se 
of the electron charge e. They further assumed [see their 
Eq. (3)] that every ion interacts with its nearest neigh- 
bors via the exponential interaction ((/)(r)/e)(^e)^, where 
(j){r) is given by Eq. ([5]). Such a semi-mean- field approxi- 
mation makes the total repulsive energy of excess anions 
larger than in our description of Sec.|lTl since the distance 
between interacting ions is smaller and this changes the 
exponential factor of (f>{r). Therefore, it seems reason- 
able that the pore capacitance evaluated in Ref. [H is 
three times smaller than in our Fig. [5l 



IV. CAPACITANCE OF A POROUS 
IMPERFECT METAL 



Thus far we have calculated the capacitance in situ- 
ations where the electrode can be considered a perfect 
metal, or in other words where the electrode has a vanish- 
ing electronic screening radius. In this section we exam- 
ine what happens to the capacitance when the electrode 
is not a perfect metal, but instead has a finite screening 
linear radius r^, given by 



e?v[^jLF) ' 



(21) 



where v{^f) is the electron density of states at the Fermi 
energy ^ip oi the electrode. 

Below we consider separately two limiting cases for : 
(i) where ^ 6, the typical thickness of the wall sep- 
arating adjacent pores, so that adjacent pores can be 
considered non-interacting, and (ii) where > h, so that 
adsorbed ions interact three-dimensionally. 

In the case where Vg <IC h, there is no interaction be- 
tween adjacent pores and the volumetric capacitance can 
still be calculated as in the previous section, by consider- 
ing the capacitance of a single pore. In this case the effect 
of finite screening radius is to shift the position of the re- 
flection plane for image charges by a distance Ts beyond 
the surface of the pore wall [22]. This reflection plane 
coincides with the "electrostatic surface" of the pore, at 
which the center of gravity of the surface charge is effec- 
tively located. That is, a charge in the center of the pore 
becomes separated from its image charge by a distance 
d + 2rs , as shown in Fig. [6] In this way the interaction 
between neighboring ions is stronger than what is given 
by Eq. ([S]) and the capacitance of the pore is reduced as 
compared to the results in Sec. |lTl One can easily calcu- 
late the effect this has on the capacitance by replacing d 
with dcs — d + 2rs in Eqs. (HI - p^ . In other words, al- 
lowing for flnite screening radius 6 in the electrodes 
has the same effect as increasing the pore width (which 
is examined in Fig. [3|) . For example, a pore with width 
d ~ 1 nm and — 0.25 nm would have d^s/d » 1.5, and 
the capacitance would correspond to the middle (blue) 
curve in Fig. [3] 











i 






FIG. 6. A schematic portrayal of the effect of finite screening 
radius ^ 6 in a conducting electrode (solid gray area). 
The reflection plane for image charges (black and white circles 
with dashed outhnes) of an adsorbed cation (black circle) is 
shifted by a distance Ts from the wall of the pore (hatched 
area). Compare to Fig. [2jb). 



In the opposite limit, > &, the wall of the pore does 
not completely screen the charge of an adsorbed cation 
and ions in adjacent pores interact with each other. In 
this limit ions interact three-dimensionally via a Yukawa- 
like potential 



(r) 



exp 



r 



(22) 



where Rs is the three-dimensional (3D) screening radius, 
determined from the volume-averaged density of states 
iy{HF)b/ib + d): 



Rs 



eoe 



d 



rs 



d 



(23) 



For d — b {a.s in Fig. [T]), we get Rs = v^r^. 

In the limit of Rs » N~^/^, where N is the three- 
dimensional concentration of ions inside the electrode, 
the electric potential $ is uniform throughout the volume 
of the electrode and is given by 



^ = eN ip{r)ATrr'^ dr 



eNRl 

£0£ 



This gives for the volumetric capacitance 



C 



(24) 



(25) 



This result was flrst derived as the volumetric capaci- 
tance of charged DNA condensates with cationic poly- 
electrolytes in salty water [23l - [25j . 

At smaller Rs < N~^/^, the discreteness of the ions 
plays an important role. In their lowest energy state, the 
ions form a correlated, 3D liquid in which they maximize 
their separation from each other while neutralizing the 
bulk charge of the electrode, as shown in Fig. [71 As in 
the previous section, we can calculate the capacitance 
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by postulating a crystalline arrangement of the ions (a 
3D Wigner crystal) and calculating the total electrostatic 
energy U. The capacitance can then be found using the 
3D analogue of Eqs. © and ([3]), namely 



V = 



dU 



1 dU 
'^dN' 



C = e^W 



[dN-^ J 



(26) 
(27) 



where V is the electrode volume. 




FIG. 7. A schematic depiction of anions (black circles) ar- 
ranging themselves within a porous supercapacitor electrode 
(solid gray area) in the limit > b. Here, ions interact 
three-dimensionally, leading to a 3D Wigner-crystal-like ar- 
rangement of ions in the ground state. 

This approach allows one to calculate C as a function 
of ion density N and as a function of voltage V. Using 
a numeric evaluation of the total energy U gives a ca- 
pacitance that can be accurately fitted to the following 
power-law form at {V - Vt)/(e/47reo£i?s) > 0.01: 



C{V) = 



1-^0.22 



V-Vt 
e/AneoeRs 



-0.58 



This expression is plotted in Fig. [5] 



c 



0.5 1 1.5 

(V-V)/(e/4jl£ eR ) 

t ./ 



FIG. 8. The volumetric capacitance C, in units of eo£/Rs, 
for the case where ionic charges obey the interaction law of 
Eq. ([22}, plotted as a function of dimensionless voltage. The 
curve is truncated at Cmax as given by Eq. (p8)l . 



Eq. ((28|) suggests that the capacitance diverges at small 
V—Vt , as in the case of 2D pores. To understand why this 
is the case, we can consider the limit where ions are suf- 
ficiently sparse that their separation is much larger than 



Rg. In this limit only the interactions between nearest- 
neighbors of the 3D Wigner crystal contribute signifi- 
cantly to the total energy. If we imagine a cubic lattice 
arrangement of ions, then we arrive at an energy 



C/~ 3V- 



.2^4/3 



■ exp 



1 



By Eqs. 



and (P7)) . this gives a voltage 

1 



V-Vt 



AireneR 



■ exp 



and a capacitance 

C{N) ~ 127reo£ViV^/3^f exp 



m/3R, 



1 



m/3R, 



(29) 



(30) 



(31) 



Combining these two relations gives an expression for the 
volumetric capacitance as a function of voltage, applica- 
ble at very small V — Vt: 



CiV) 



R^'^V-Vt 



In" 



e/ineoeRs 
V-Vt 



(32) 



Eq. ((32|) implies that at small V — Vt, where ions are 
sparse, the capacitance can be much larger than the 
mean- field result of Eq. (f25|) . This growth in the ca- 
pacitance is driven by the vanishing interaction [Eq. (j22p ] 
between discrete, correlated ions. The maximum value of 
the capacitance occurs at (V—Vt) — > and is determined 
by thermal effects. This maximum can be estimated, as 
in the previous section, by making a virial expansion of 
the free energy 



F~F,d+ YN^kBTB{T) - Q{V - V] 



(33) 



The value of the capacitance in this limit, as in Eqs. (jl7|) - 
(IT51) . is inversely related to the virial coefficient B{T). At 
not too small i?s, such that Rg > a/2, 



2T:eoeRs 
fB{T) 



(34) 



Here, T is a dimensionless temperature normalized to the 
interaction between two charges at a distance Rg : 



kBT 



d 



e^/ATTEaeR, 
The virial coefficient B{T) is calculated as 

oo 

e0(r) 



B{T) (^1 - exp 





2'kRI 



47rr dr 



l + ln[l/T] + -ln^[l/T] 



SO that Eq. (IMl) can be written 



£oe 



T(l + ln[l/T] + iln3[l/r] 



i?2 



(35) 

(36) 
(37) 

(38) 
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If one wishes to formulate an approximate prediction 
for the capacitance at arbitrary values of screening ra- 
dius, including values of rg that are comparable to b, then 
one may evaluate separately the capacitance based on ap- 
proach (i), where the 2D pore thickness is renormalized 
as d "> doff, and approach (ii), where ions interact three- 
dimensionally, and then take the smaller value. Since 
intra- and inter- pore interactions contribute additively 
to the total energy, these create series contributions to 
the capacitance, so that as a zero-order approximation 
one can take the smaller of the two capacitances. The 
result of this process is shown in Fig. ^ which consti- 
tutes a prediction for the volumetric capacitance C at 
arbitrary voltage and electrode screening radius. 
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FIG. 9. (Color online). A contour plot of the volumetric ca- 
pacitance C as a function of voltage V and electrode screen- 
ing radius for charges in an electrode with planar pores of 
width d separated by walls of thickness fe = d at room tem- 
perature. Contours are labeled by their values of C in units 
of eoe/d? . The dashed line separates the regions of validity 
of the two theories presented in this section: (i) where pores 
are non-interacting and the pore width can be renormalized 
according to d — >■ — d + 2rs (lower region), and (ii) where 
charges interact three dimensionally with the interaction law 
of Eq. (|22l) (upper region). For e = 2 and d = 1 nm, the unit 
of volumetric capacitance eos/d^ = 18 F/cm^ and the unit of 
voltage e/4neosd = 0.7 Volts. 



V. NONLINEAR SCREENING IN GRAPHITE 
(CARBON) SUPERCAPACITORS 

The preceding sections outline a general theory for the 
volumetric capacitance of a supercapacitor made with 
electrodes that can either be considered metallic or can 
be characterized by some linear screening radius rg. In 
this section we discuss specifically the case of graphite 
electrodes, which is among the most comm only -studied 
materials for supercapacitor devices 0, H, 0, ES] • 

In graphite, the Fermi level density of states v{hf) is 
actually relatively small, so that the linear screening ra- 
dius Ts w 0.8 nm cannot be considered much smaller than 
the spacing between pores [26j . For example, if d = 6 = 1 



nm, then Fig.|9]would seem to imply a capacitance on the 
order of eoe/(P ~ 18 F/cm'^. Experiments with graphite 
electrodes, however jvield a capacitance five times larger 
than this value [1, 0], suggesting that graphite screens 
over a much smaller distance than and effectively be- 
haves as a good metal. 

This apparent discrepancy can be resolved if one recalls 
that the density of states in graphite is close to i^^fJ^p) 
only in a narrow range of energies, beyond which it in- 
creases linearly with energy on both sides of the Fermi 
level dl] . Such variation of the density of states suggests 
that screening by graphite is non-linear even when a rel- 
atively small electric field is applied to the surface of the 
pore, and that therefore the screening properties of the 
electrode material cannot be characterized by a constant 
linear screening radius r^. 

In order to estimate the distance over which the ions' 
potential is screened, one can consider the problem of a 
uniform applied electric field Eq orthogonal to the basal 
plane of graphite. It has been shown [23] that in this case 
the magnitude of the electric field E{z) decays with the 
distance z beyond the graphite interface as 



where 



Eo 



(1 + z/2zor 



ATTEQea^c^EQ 



1/3 



(39) 



(40) 



is the centroid location of the counter-charge in the 
graphite (the non-linear screening radius), c « 0.34 nm is 
the distance between graphite planes (graphene sheets), 
and a « 2.2/e is the effective fine structure constant of 
graphene. (See also recent discussions of screening in 
graphene multilayers in Refs. and[29h. 

The implications of this result for the problem of 
screening of adsorbed ions in a graphite pore can be seen 
as follows. When the area density of ions inside a pore 
is large enough that na^ ~ 1, these ions can be said to 
produce a roughly uniform electric field at the wall of the 
pore whose strength is Eg ~ ej^e^ea^ . Inserting this re- 
lation into Eq. pO|) yields a screening distance zq ~ 0.75c, 
which suggests that the field is entirely screened within 
the first graphene layer. In other words, at dense ion 
filling the electric field between ions in the pore does not 
penetrate beyond the first graphene layer. Therefore, de- 
spite its relatively low density of states v(jxp), graphite 
may be treated as a metal at not-too-small ion densities 
nd^ > 0.5 (note that no} > 0.5 occupies the majority of 
the voltage range in Fig. [3]). 

We can also discuss what happens with the volumetric 
capacitance when the density of ions is smaller (V — Vt 
decreases), so that the electric field Eq produced by the 
ions becomes weaker. When V — Vt is made moderately 
small, the effective density of states decreases, the non- 
linear screening radius zq grows, and the nonlinear ca- 
pacitance decreases. One can show using Eq. that 
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the capacitance C oc {V — Vt)^^^ . Eventually, at small 
V—Vt, the effective density of states saturates at the level 
of v{iif), so that the screening radius becomes constant 
and equal to the linear screening radius ~ 0.8 nm. In 
this limit, the volumetric capacitance is relatively small 
and is given by the theory of Sec. IIV) as was already 
discussed in the beginning of this section. 

VI. CAPACITANCE OF A CRYSTALLINE 
ASSEMBLY OF METALLIC SPHERES 

So far we have restricted our discussion to the elec- 
trode geometry shown in Fig. [TJ In practice, such elec- 
trodes with parallel planar pores are difficult to make. In 
many cases supercapacitor electrodes are simply a ran- 
dom assembly of conducting particles, arranged so that 
the particles form an infinite, conducting cluster through 
which electrons can percolate while the pores in this clus- 
ter form a separate percolating space through which the 
ionic liquid can freely pass. For such cases the model of 
Fig. [I] is a strong idealization. In this section we would 
like to briefly discuss another idealized electrode struc- 
ture which captures the three-dimensional character of 
pores. 

Consider an assembly of metallic nanospheres, each 
with the same radius i?, arranged so that they form a 
cubic lattice with nearest neighbor spheres touching each 
other. As in previous sections, we imagine that this crys- 
talline film is deposited onto a metallic contact plate, 
connected to a voltage source, and immersed in an ionic 
liquid. A voltage V is applied between the contact plate 
and the bulk of the ionic liquid. If the diameter of the 
ions within the ionic liquid is small enough, then this 
arrangement produces an effective supercapacitor elec- 
trode, where ions may percolate through the spaces be- 
tween conducting spheres and neutralize the electronic 
charge provided by the voltage source. 

In order to calculate the capacitance of this electrode. 



we first analyze the interaction between two ions that en- 
ter into the bulk of the electrode. This can be done by 
calculating the potential as a function of distance pro- 
duced by a single ion in the center of a pore deep inside 
the electrode bulk. We calculate this potential numeri- 
cally using the relaxation method for solving the Laplace 
equation, where each of the conducting spheres is held at 
zero potential. We find that potential decays exponen- 
tially with radial distance from the ion with a character- 
istic screeening length Rg = Rx (0.26±0.01). The reason 
for this sharp decay is the same as for the decay of the 
potential in 2D, slit-like pores [Eq. ©]: electric field lines 
emanating from the ion are adsorbed by the surface of 
nearby conducting spheres, and the number of these field 
lines that survive by passing through the narrow, tor- 
tuous pores between spheres decays exponentially with 
distance. 

Once the interaction law is known, one can calculate 
the capacitance in a way similar to the analysis of Sec. 
IIVI We arrive then at a relation C{V) which has the 
maximum given by Eq. (I38|) . For i? = 4 nm, T = 300 
K, and e — 2 this relation gives a volumetric capacitance 
C^^^ = l.lsoe/Rl = 25eoe/R'^ « 28 F/cm^. Remark- 
ably, in this arrangement the capacitance per sphere is 
roughly 16 times larger than the capacitance of a single, 
isolated sphere in a medium with dielectric constant e. 

One can reach even larger volumetric capacitance if 
the spheres are densely packed rather than arranged in a 
cubic lattice. Reducing the radius R of the spheres also 
sharply increases the capacitance. In our next publica- 
tion we will explore these mechanisms for increasing the 
capacitance by combining MC modeling with the analyt- 
ical estimates presented in Sec. IIVI 
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